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A fully self-consistent model of the neutron star inner crust based upon models of the nucleonic 
equation of state at zero temperature is constructed. The results nearly match those of previous 
calculations of the inner crust given the same input equation of state. The extent to which the 
uncertainties in the symmetry energy, the compressibility, and the equation of state of low-density 
neutron matter affect the composition of the crust are examined. The composition and pressure of 
the crust is sensitive to the description of low-density neutron matter and the nuclear symmetry 
energy, and the latter dependence is non-monotonic, giving larger nuclei for moderate symmetry 
energies and smaller nuclei for more extreme symmetry energies. Future nuclear experiments may 
help constrain the crust and future astrophysical observations may constrain the nuclear physics 
input. 
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I. INTRODUCTION 

The inner crust of a cold neutron star can be defined 
as the region between the density where neutrons drip 
out of nuclei (about 4 xlO 11 g/cm 3 ) and the density for 
the transition to homogeneous nucleonic matter at about 
half of the nuclear saturation density. This region is sen- 
sitive to the nuclear physics input because the nature of 
the crust is determined by the structure of neutron-rich 
nuclei and the energetics of the surrounding dripped neu- 
trons. In this work, the dependence of the description of 
the neutron-rich nuclei and the dripped neutrons on the 
equation of state (EOS) of homogenous nucleonic matter 
is examined. 

The inner crust is of broad interest because a large va- 
riety of astrophysical observations are dependent on and 
sensitive to the properties of the neutron star crust. One 
recent motivation is the suggestion that the giant flares in 
Soft Gamma-Ray repeaters trigger seismic events in the 
neutron star crust and are sensitive to the shear modulus 
of the crust crust [H, 0, The shear modulus, in turn, 
is sensitive to the composition of the neutron star crust 
and the relative magnitude of the proton and neutron 
numbers of the nuclei in the inner crust. The moment of 
inertia of rotating neutron stars is also sensitive to the in- 
ner crust and depends on the transition density between 
the crust and the core Q. Neutrino and photon opac- 
ities are also sensitive to the properties of the nuclei in 
the inner crust. For example, neutrino-nucleus scatter- 
ing, which scales like A 2 , is the most important neutrino 
process during the lepton-trapped phase of a Type II su- 
pernova (see Ref. Q for a recent review). Finally, the 
cooling and evolution of neutron star crusts depends on 
the both the size of the crust Q and by its transport 
properties @, [1[ , which are both related to the composi- 
tion. These astrophysical connections motivate the study 
of the magnitude of the uncertainty of the properties of 
the inner crust which come from present uncertainties in 
the nuclear physics inputs. 

In this article, several models of the neutron star inner 



crust are constructed systematically using inputs from 
with the current experimental information while allowing 
the range of uncertainty allowed due to the uncertainty 
in the EOS of homogeneous matter. Of particular impor- 
tance, is that the symmetry energy is varied in both the 
description of the nuclei and the description of the neu- 
tron matter at the same time. The composition depends 
on the symmetry energy, but is nearly independent of the 
compressibility This means that astrophysical observa- 
tions which arc connected to properties in the crust can 
constrain the nuclear symmetry energy. 

The inner crust is quite sensitive to the EOS of neutron 
matter at sub-saturation densities. At sufficiently low 
densities, neutron matter is somewhat well understood 
because three-body interactions are small, and the two- 
body neutron-neutron interaction is strongly constrained 
by the experimentally measured neutron-neutron scat- 
tering phase shifts @ (see also the review in Ref. [Iol|). 
Many of the currently available EOSs, however, do not 
respect this understanding of low-density neutron matter 
because they are fit to the properties of nuclei which are 
more sensitive to matter near saturation densities. The 
neutron matter EOSs used in this work are designed to 
have a realistic behavior below the saturation density, 
within the precision required for the description of the 
crust. 



II. THE MASS MODELS 

While microscopically-based models of the nuclei are of 
great interest because they can disentagle important ef- 
fects which are not easily treated in a classical approach, 
a microscopic approach can also make it more difficult to 
understand the physical principles which guide the na- 
ture of the inner neutron star crust. In addition, it is not 
clear that a classical approach is significantly less effective 
at estimating the magnitude of uncertainties originating 
in the nuclear physics input (it may even be more effec- 
tive). In any case, since the purpose is only to estimate 
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the uncertainties from the nuclear physics input to the 
EOS, a liquid-drop model quite similar to that described 
in Refs. [ill, EH is used. More microscopic models for the 
crust have been developed (see the pioneering work of 
Ref. [IH and recent efforts in Refs. [lj, [Hj]) and it is ex- 
pected that these results on the sensitivity to the EOS of 
homogeneous nucleonic matter will apply to some extent 
in these models as well. 

The liquid-drop model for this work consists of a bulk 
energy contribution which is determined from the EOS 
of homogenous nucleonic matter together with surface 
and Coulomb contributions. This will be compared to 
the finite-range droplet model described in Ref. [l6| and 
used in Ref. [l7[ . 

The binding energy per baryon of a nucleus with pro- 
ton number Z and atomic number A is given by 

B(Z, A) /A = B hulk (n n ,n p )/A + aB(n n ,n p ) l^j) 

+C ECoulomb/?'* (1) 

where n n and n p are the average neutron and proton 
densities inside the nucleus with the given Z and A. The 
binding energy of bulk matter, -Bbuik (about —16 MeV in 
isospin symmetric matter) is given by 

-Bbuik = — [ £ ( n n, n p ) - n n m n - n p m p ] (2) 

where m n and m p are the neutron and proton masses 
(which is taken to be 939 MeV) , n = n n + n p is the aver- 
age baryon number density in the nucleus, and e{n ni n p ) 
is the energy density of homogeneous matter evaluated 
at the given neutron and proton density. The expression 
for e may be given by any EOS of homogeneous matter 
and several different models are employed. Note that the 
energy of the dripped neutrons which is added later will 
always be determined with the same EOS as is used to 
describe the bulk part of the nuclear energy 

The average baryon density will be determined from 

n = n n + n p = n + nil (3) 

where 1 = 1 — 2Z/A. The parameter no is analogous to 
the saturation density of nuclear matter and is expected 
to be near 0.16 fm -3 . The parameter n\ subsumes (in 
a very schematic way) two effects: the decrease in the 
saturation density with the isospin asymmetry and the 
increase in the saturation density due to the Coulomb in- 
teraction. These effects are both explicitly present in the 
finite range droplet model (see Eq. 49 of Ref. [13 )■ The 
decrease in the saturation density with isospin asymme- 
try is typically larger and thus n\ is always negative in 
these models. 

The individual average neutron and proton number 
densities are given by 

n n = n(l + 5)/2 

n p = n(l-<5)/2 (4) 



and the density asymmetry S = l — 2n p /(n n + n p ) is given 
by 5 = (I where £ is a constant parameter of the model. 
Neutron and proton radii ( "squared-off " radii, not root- 
mean-square radii) are given simply by 47m„i? 3 = 3N, 
and 4:Trn p R p = 3Z. The presence of a neutron skin is 
determined from £. If £ is unity, then all nuclei have no 
neutron skin (R n = R p ), while if £ is less than unity, 
then all nuclei with N > Z will have a neutron skin 
(R n > R p ). 

The surface energy contribution is proportional to the 
surface tension a, A 2 / 3 (the surface energy scales as A 2 / 3 
so that the surface energy per baryon scales like A -1 / 3 
as in Eq. [I]), and a unitless function B. Typically this 
latter function is quadratic in the isospin symmetry 

B(n n ,n p ) = 1 - a s S 2 (5) 

where as is a positive parameter representing the surface 
symmetry energy. This is essentially the approach taken 
in Ref. [16| . For a neutron star inner crust model, this can 
be modified to ensure that the surface energy vanishes in 
the limit 8 — > 1 as it must. One possible approach (and 
the one used here) is that from Ref. [ljj 

ni \ 16 + b (r\ 

BK " ?) = [iFWiF^ () 

where x = n p /n and b is a simple function of the pa- 
rameter as and is related through as = 96cr/(6 + 16). 
This is an approximate scheme for taking into account 
the isospin properties of the surface energy which may 
suffice for the present purpose, but note the more the 
more detailed discussion in Ref. [H[ . In particular there 
is still an unresolved ambiguity associated with how the 
surface energy is handled as discussed in this reference. 
The slope of the correlation between the surface symme- 
try energy, as and the symmetry energy at the saturation 
density depends on the mass formula used. This model, 
like all other present models for the neutron star crust, 
effectively chooses a particular slope for this correlation. 

The Coulomb energy density of a 3-dimensional 
droplet of protons can be written [ll|, [l9[ (modulo an 
overall factor of x, the volume fraction of matter present 
in nuclei, which is included later), 

^Coulomb = Y n l e2R l ( 2 - 3 * 1/3 + X) (7) 

where e 2 is the usual Coulomb coupling ~ he/ '137. In the 
final term in parenthesis, the first term corresponds to the 
standard Coulomb contribution, the second term corre- 
sponds to the "lattice contribution" [2(j in the Wigner- 
Seitz approximation, and the last term to a further finite- 
size correction relevant at higher densities when x is com- 
parable to unity. Note that this last term is quite impor- 
tant near the crust-core transition and tends to delay 
the transition to nuclear matter to higher densities. The 
Coulomb contribution is multiplied by a parameter C to 
take into account the fact that the proton density does 
not fall off sharply at a finite radius and this surface 
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diffusiveness maybe dependent on the input symmetry 
energy. This parameter will always be nearly unity. As 
noted in Ref. [2l|, the Wigner-Seitz approximation fails 
when describing the low-temperature transport proper- 
ties, but will suffice for describing the composition and 
the equation of state as done here. 

In summary, there are six free parameters in this model 
(outside of the input equation of state of bulk nuclear 
matter, which is a kind of parameter in itself) are the 
surface tension in MeV/fm 2 , a, the surface symmetry 
energy as, the correction factor to the Coulomb energy, 
C, the asymmetry parameter £, and the central density 
parameters, hq and ri\ which are expressed in units of 
fm . These six parameters will be fit to experimental 
masses for each input EOS. 



III. THE EQUATIONS OF STATE 

The EOS from Ref. 01 (APR) is used, which was ob- 
tained from variational chain summation calculations of 
the equation of state using a realistic nucleon-nucleon in- 
teraction. Also, a "typical" relativstic field-theoretical 
model is utilized (review in Ref. [H[), NL4 [24j which 
was fit to nuclei. In order to compare with the model of 
Ref. [H the Skyrme [p model SLy4 ^ is usecL and 
in order to compare with the model from Ref. [28[ the 
Skyrme model SkM* [H is used. 

APR is expected to be particularly good for neu- 
tron matter at low densities, because it is directly com- 
puted from an interaction which reproduces the two-body 
nucleon-nucleon phase shifts. The model SLy4 also has a 
good neutron matter EOS because it was fit to both nu- 
clei and low-density neutron matter. The NL4 and SkM* 
models were only fit to nuclei and low-density neutron 
matter are less constrained. SkM* happens to have a 
neutron matter EOS which is somewhat closer to APR 
than NL4. Like the SLy4 interaction, rclativistic mod- 
els are also able to reproduce, at some level, the more 
accurate low-density neutron matter EOS found in APR 
and SLy4, as was demonstrated by the RAPR model in 
Ref. and the FSUGold model 

In order to examine the importance of having an accu- 
rate EOS for low-density neutron matter, the low-density 
neutron matter EOS of NL4 is modified and compared to 
the original. Three new models are constructed. The first 
model, NL4Q, is a modification of NL4 which treats the 
symmetry energy at low densities to be exactly quadratic. 
This approximation is quite good at lower densities, and 
the NL4 crust is nearly indistinguishable from the NL4Q 
crust. The other models, NL4QN and NL4QN2, are ver- 
sions of NL4Q which reproduces the neutron matter EOS 
of APR at densities below a specified density, 
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1 -|- e {n t -n)/v 



(8) 



where n t is 0.08 (0.04) fm" 3 and v is 0.0105 (0.016) 
fm" 3 for model NL4QN (NL4QN2). Both NL4QN and 
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TABLE I: The values of 771 and 772 (c.f. Eq. |9]), parame- 
ters controlling the neutron matter equation of state in the 
schematic models and the corresponding values of the sym- 
metry energy at saturation density and the exponent 7. In 
the schematic models, 770 is fixed at 0.5. 



NL4QN2 have pressures and neutron chemical potentials 
which monotonically increase with density. The results 
below show that these two models give signficantly dif- 
ferent results for the inner crust. 

Finally, several schematic models of the EOS are con- 
structed so that the effect of the compressibility and the 
symmetry energy can be examined. The schematic EOS 
for neutron matter is 
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(9) 



where kp = (37T 2 ?!) 1 / 3 is the neutron Fermi momentum 
and the last term is just the free Fermi gas energy den- 
sity. The first two terms inside the square brackets are 
designed to reproduce the expectation from equations of 
state at low densities obtained from two-body potentials 
which reproduce the experimental phase-shift data on 
neutron-neutron scattering 0, [32j. It is expected that 
the interacting neutron matter EOS is about half the free 
Fermi gas energy at kp^ n = 0.5 fm -1 . The term propor- 
tional to k F A n , qualitatively reproduces this low-density 
behavior and the other parameters 771 and 772 can be ad- 
justed for densities near the saturation density where the 
EOS is more uncertain. 

It is useful to connect this description with the more 
traditional description neutron matter in terms of a 
symmetry energy with the form E sym = A(n/no) 2 / 3 + 
B(n/no) 7 . For an effective mass of about 0.7 M, A is 
about 17 MeV, and then B and 7 dictate the magnitude 
of the symmetry energy at the saturation density and the 
density dependence of the symmetry energy, respectively. 
The base model schematic model, "Sch" , has values of rji 
and 772 appropriate for B = 14 MeV (a symmetry energy 
of 31 MeV) and 7 = 0.9. These values arc given in Ta- 
ble I, as well as the values of 771 and 772 for the schematic 
models whose symmetry energy is different from the base- 
line model. The variation of the value of the symmetry 
energy at saturation density between 28 and 34 MeV is 
consistent with the observation that most modern equa- 
tions of state fall within this range. When the symmetry 
energy is taken to be a pure power law (the A=0 limit), 
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the limit of the variation of 7 (0.6 < 7 < 1.1) is inferred 
from the experimental information from intermediate- 
energy heavy- ion collisions [HI, [HI]. The models Schyl 
and SCI172 are constructed by fitting a symmetry energy 
with the given exponent with A=0 to express the ex- 
pected range. The extraction of values of 7 from heavy- 
ion collisions is non-trivial, and there may be systematic 
uncertainties that are not yet understood. These uncer- 
tainties would mean that the range of variation presented 
here is overly conservative, and that the true range might 
be larger. 

The schematic equation of state for nuclear matter is 



E nuc /A = AI + B 
K' 



K ( ,2 
— (n - n ) 



162ng 



18n 



(n - no)' 



(10) 



where M is the nucleon mass, B is the binding energy, 
no is the saturation density, K is the compressibility 
and K' is the "skewness" . Isospin asymmetric matter is 
computed assuming that the symmetry energy is exactly 
quadratic in the isospin asymmetry S (thi s approxima- 
tion may fail at high density, see Ref. [35|). Note that 
varying the compressibility in this model is not precisely 
equal to varying the quantity which might be obtained 
from giant resonances, as the latter are only sensitive to 
the equation of state in the neighborhood of saturation 
density whereas the compressibility is applied to nuclear 
matter at all densities below the saturation density, thus 
the variation in the compressibility is a bit larger than 
that recently suggested in Refs. [3g, H3, [Hj]. The base- 
line schmatic model "Sch", has a binding energy of —16 
MeV, a saturation density of 0.16 fm -3 , a symmetry en- 
ergy of 31 MeV, a compressibility of 230 MeV. In all of 
the models, the skewness parameter is fixed by ensuring 
that the energy per baryon of nuclear matter vanishes 
at zero density, as it ought. In addition to variations 
of the symmetry energy as discussed above, Two mod- 
els "SchK210" and "SchK250", are constructed to be the 
same as the baseline model, except that they have differ- 
ent compressibilities. 

A survey of some of the equations of state is given in 
Fig. [TJ The upper right panel shows how neutron matter 
in the NL4QN and NL4QN2 models interpolates between 
APR at low densities and the normal NL4Q model at 
higher densities, while leaving the nuclear matter EOS 
unmodified. 
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FIG. 1: A survey of the equations of state used in this work. 
Plotted in each panel are the energy per baryon of nuclear 
matter (lower set of curves) and neutron matter (upper set 
of curves) as a function of baryon density, ub- The upper 
left panel shows the models APR, SLy4, and SkM*, the up- 
per right panel shows the model NL4Q (based on NL4) and 
shows how neutron matter was modified to match APR at 
low densities, and the lower panels show how the symmetry 
energy was modified in the schematic equation of state. 
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FABLE II: The nuclear mass fits corresponding to the models 
described in the text. The values of a are given in MeV/fm 2 
and no and ri\ are given in fm -3 . 



IV. THE MASS FITS 

The liquid drop model is fit to the experimental nuclear 
masses from Ref. [39| using the fitting formula 
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(11) 



where TV is the number of nuclei, and M exp and M th 
are the experimental and theoretical values of the mass 
excess. Ref. [17[ points out that this fitting formula can 
be improved and that it overestimates the actual model 
error, but these considerations will not be important at 
the level of the results presented here. The fitting results 
are given in Table II. The fitting results for NL4QN and 
NL4QN2 are not given because they were found to be 
nearly equal to those from NL4Q. 
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The mass fit is performed by minimizing Arms f° r 
all the experimentally measured mass excesses from 
Ref. For each nucleus, this involves computing the 
neutron and proton densities using Eqs. [3] and 31 com- 
puting the bulk energy from the EOS of homoegenous 
matter at these densities, then inserting this bulk energy 
into the nuclear mass formula to compute the mass excess 
using Eq. [TJ 

As expected, there is a correlation between the sur- 
face symmetry energy and the symmetry energy as shown 
by the increase in a$ when going from model SchS28 to 
model SchS34. This is also the reason why NL4 gives a 
larger value of as than the other models. The value of ag 
is nearly unchanged by modifying 7, which changes the 
density dependence of the symmetry. The values of no in 
Table [TV] are not quite equal to the saturation density for 
homogeneous nuclear matter, and this can be attributed 
to finite-size effects not captured in Eq. [3j The other 
parameters are nearly unchanged between models except 
for fii which is also sensitive to the symmetry energy, as 
well as the Coulomb interaction. 



V. THE CRUSTS 

In order to determine the composition and properties 
of the crust, the energy at a fixed density as a function 
of the proton number and atomic number of nuclei, and 
the number density of dripped neutrons, n rh drip, is mini- 
mized. The energy of matter in the neutron star crust is 
given by 

e{Z, A, n„ idrip ) = (n„ + n p )xB(Z,A)/A + 

(1 - X^dripOVdrip) + £cl(™ e )(12) 

This energy is minimized over the three parameters Z, A, 
and n n! drip at each density. The volume fraction of mat- 
ter inside nuclei, x, is determined from the relation 



tib = x(n„ + n p ) + n„ idri p(l - \) 



(13) 



The inner crust implied by the models NL4Q, NL4QN 
and NL4QN2 are compared in Fig. [3] While the actual 
number density of dripped neutrons is not strongly mod- 
ified by modifying the neutron matter EOS, the nuclear 
size is modified by 50 % or more. The larger energy cost 
of creating neutron matter in with a more realistic neu- 
tron matter EOS is reflected in moving neutrons into nu- 
clei so as not to pay the energy cost. Because the largest 
difference is between the models NL4Q and NL4QN2, it 
is clear that most of the dependence on the low-density 
neutron matter EOS lies at densities below 0.04 fm~ 3 . 

In order to show the effect due to changing the mass 
model, Fig. [3J shows the results for the APR EOS with 
the two different models, the liquid drop model and the 
FRDM. The results are qualitatively the same but quan- 
titatively different. The form of the mass model remains 
a significant uncertainty in the nature of the neutron star 
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FIG. 2: A comparison of the composition of the neutron star 
crusts in models NL4Q, NL4QN, and NL4QN2, designed to 
demonstrate the sensitivity of the composition to the equation 
of state of low-density neutron matter. The proton nuber, Z, 
the atomic number A (left axis) and the number density of 
the dripped neutrons n ni dri P are displayed (right axis). 
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FIG. 3: A comparison the composition given the two different 
mass models. The curve labeled LDM (solid lines) is obtained 
from Eq. [1] and the curve labeled FRDM (dashed lines) is 
obtained from Ref. [13] ■ The left axis is for A and Z and the 
right axis is for n n ,drip- 
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FIG. 4: A comparison of the present work (solid lines) with 
the results of Ref. [25l ] (dashed lines) using the same input 
equation of state, SLy4. The left axis is for A and Z and the 
right axis is for n n: drip- 



crust, and is comparable to the other uncertainty ob- 
tained from the symmetry energy as described below. 

Comparisons of the present model to that of Rcfs. (2f| 
[H[ are given in Figs. 2] and [51 The same input EOS for 
homogeneous matter as the original reference is used in 
both cases. The results agree qualitatively with the afore- 
mentioned works. The remaining differences lie within 
the nuclear mass formula used, and they are within the 
range of variation which is suggested by Fig. [3] This 
model (like all other models of the neutron star crust 
presently available) cannot precisely predict the compo- 
sition of the inner crust. Nevertheless, it is qualitatively 
correct and is thus useful for estimating the uncertainties 
due to the input EOS of homogeneous matter. 

To compare the effect of the uncertainty in the sym- 
metry energy, Fig. [5] shows the composition for the neu- 
tron star crust as a function of density for the schematic 
equations of state with different symmetry energies. 
The naive expectation is that a stronger symmetry en- 
ergy tends to encourage nuclei to become more isospin- 
symmetric. This is coupled, however, with the fact that 
an increased symmetry energy will also raise the energy 
cost for the dripped neutrons. These two effects together 
could force larger, more symmetric nuclei, but this also 
affects the Coulomb and surface energy contributions. 
The variation of the composition with the value of the 
symmetry energy is not so clear, as the baseline model 
predicts larger nuclei than either models with smaller or 
larger values of the symmetry energy. 

In order to disentangle this result, more detailed results 
for schematic models with different symmetry energies 



FIG. 5: A comparison of the present work (solid lines) with 
the results from Ref. [H (labelled "JML" and plotted with 
dashed lines) using the same input equation of state, SkA/* . 
The sharp dropoff in JML at high- densities is due to the tran- 
sition to homogeneous nucleonic matter. The slightly jagged 
nature of the JML results is due to the naive interpolation 
employed in this work and is not necessarily present in the 
original table. The left axis is for A and Z and the right axis 
is for 7i nidrip . 
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FIG. 6: A comparison of the composition of the crust given 
different symmetry energies. The bold solid line is the baseline 
model, the dashed and dotted lines show the variation in 7, 
and the dashed-dotted and thin solid line give the variation 
in the magnitude of the symmetry energy at the saturation 
density. 
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FIG. 7: A comparison of the composition of the crust given 
different values of the symmetry energy at the saturation den- 
sity at a fixed density of ub = 0.01 fm -3 . The top panel 
gives the total binding energy per baryon, and the separate 
contributions from dripped neutrons ("n"), electrons ("e"), 
and nuclei ("Nuc"). The bottom panel shows the neutron 
and proton number of nuclei as well as the volume fraction, 
X- 

are given in Fig. [7] at a fixed density of ub = 0.01 fm -3 . 
Beginning with the larger symmetry energy (with a value 
at saturation of 34 MeV) and proceeding downward, the 
expected result is obtained: lower symmetry energies al- 
low the system to create more isospin-asymmetric nuclei. 
At low enough symmetry energies, however, this becomes 
too costly as the electron contribution to the energy in- 
creases (the proton number decreases, but the volume 
fraction occupied by nuclei increases, thus the electron 
density must increase). Instead, the system reponds by 
moving neutrons out of the nuclei, which lowers the elec- 
tron contribution, even though it increases the contri- 
butions from nuclei and the dripped neutrons. This is 
allowed, in part, because the nuclei are able to maintain 
a relatively constant energy. They can do this because 
the surface and Coulomb energy cost is cancelled by the 
bulk energy gain which results from making nuclei with 
a larger (in absolute magnitude) bulk binding energy. 

Finally, Fig.[S]summarizes the pressure as a function of 
the baryon density. The upper left panel shows the vari- 
ation from the different mass formulas, which is larger at 
lower densities. The upper right panel shows the results 
for SkM* and SLy4. The lower-left panel shows the vari- 
ation allowed by the symmetry energy, and variations of 
up to a factor of two in the pressure are implied by the 
uncertainty in the symmetry energy The pressure ap- 
pears sensitive to the magnitude of the symmetry energy 
and its dependence on density. Finally, the lower-right 



FIG. 8: A survey of the pressure of the crust as a function of 
density, scaled by the baryon density, n^ 3 . The curves are 
nearly flat where the adiabatic index is expected to be nearly 
4/3 since the baryon density nearly scales with the energy 
density. 



panel shows the pressure for the NL4Q-related models. 
Note that for model NL4Q, the anomalously small EOS 
of low-density neutron matter underestimates the pres- 
sure in the inner crust. All of the models are connected 
(sometimes discontinuously) to the EOS from Ref. [2(| 
at low densities thus giving the scatter in the pressure 
at the lowest densities given in this figure. Particularly 
interesting is that the pressure at the crust at the higher 
densities is nearly independent of the mass model (as 
shown in the upper left panel), which may indicate that 
astrophysical observables which arc sensitive to the pres- 
sure of the crust rather than the composition are good 
probes of the nuclear symmetry energy, assuming that 
the description of low-density neutron matter is correct. 

The transition density to homogeneous nuclear matter 
is computed by noting the density at which the energy 
per baryon of nuclear matter becomes smaller than that 
of the heterogeneous phase. The results for the transi- 
tion density arc given in Table III. The transition densi- 
ties from this work are given in the second column and 
third column contains the transition densities for previ- 
ous works with similar input EOSs for comparison. These 
transition densities may be underestimates, as the abil- 
ity of nuclei to deform slightly will decrease the energy 
of the heterogeneous phase and thus increase the tran- 
sition density. This will be addressed in further work. 
Note that, as in the composition discussed above, the 
transition densities depend non-trivially on the symme- 
try energy. 
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Model 


rit fm~ 3 




APR 


0.0522 




SkM* 


0.0434 


0.045 [28] 


SLy4 


0.0669 


0.076 [25] 


NL4Q 


0.0344 




NL4QN 


0.0409 




NL4QN2 


0.0333 




Sch 


0.0584 




SchK210 


0.0585 




SchK250 


0.0591 




SchS28 


0.0368 




SchS34 


0.0416 




Sch 7 l 


0.0676 




Sch 7 2 


0.0641 





TABLE III: Caption here. 



VI. CONCLUSIONS 

The composition of the neutron star crust is still par- 
tially unknown, due to uncertainties in the nuclear mass 
formula and the equation of state. The composition (and 
to a lesser extent, the overall pressure) is quite sensi- 
tive to the equation of state of low-density neutron mat- 
ter, and the nuclear symmetry energy, both its magni- 
tude and its density dependence. The dependence of the 



composition on the symmetry energy is not monotonic, 
as models with moderate symmetry energies can have 
larger nuclei than models with lower or higher symmetry 
energies. To the extent to which neutron stars depend 
on the composition, this means that it is important to 
explore the full range of variation in the crust allowed 
by the present knowledge of the input nuclear physics, 
while ensuring that the EOS is constrained by what is al- 
ready known about the EOS of low-density neutron mat- 
ter. Nuclear experiments will continue to provide better 
constraints on the symmetry energy, including from the 
PREX experiment [10, [H| to measure the neutron skin 
thickness of lead at Jefferson Lab and from intermediate- 
energy heavy-ion collisions as has been done in Ref. [H] . 

It remains to be seen if these results persist in the more 
microscopic models which include pairing, corrections be- 
yond the Wigner-Seitz approximation, long-range corre- 
lations, and better treatments of the nuclear structure. 
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